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Ground state and bias current induced rearrangement of semifiuxons in O-rr long 

Josephson junctions. 

E. Goldobing D. Koelle, and R. Kleiner 
Physikalisches Institut - Experimentalphysik II, Universitat Tubingen, 
Auf der Morgenstelle 14, D-72076 Tubingen, Germany 
(Dated: February 1, 2008) 

We investigate numerically a long Josephson junction with several phase 7r-discontinuity points. 
Such junctions are usually fabricated as a ramp between an anisotropic cuprate superconductor like 
YBa2Cu307 and an isotropic metal superconductor like Nb. From the top, they look like zigzags 
with 7r-jumps of the Josephson phase at the corners. These 7r-jumps, at certain conditions, lead to 
the formation of half-integer flux quanta, which we call semifiuxons (SF), pinned at the corners. 
We show (a) that the spontaneous formation of SFs depends on the junction length, (b) that the 
ground state without SFs can be converted to a state with SFs by applying a bias current, (c) that 
the SF configuration can be rearranged by the bias current. All these effects can be observed using 
a SQUID microscope. 

PACS numbers: 74.50.+r, 85.25.Cp 74.20.Rp 

Keywords: Long Josephson junction, sine-Gordon, half-integer flux quantum, semifluxon, 0-pi-junction 



I. INTRODUCTION 

Due to the specific order parameter symmetry of high- 
T c cuprate superconductors it is possible to fabricate and 
study so-called and tt Josephson junctions (see Ref. |l| 
and references therein). These junctions can be fabri- 
cated in various geometries and open new perspectives 
for Josephson electixinics (digital circuits, fiuxon devices, 
quantum bits, etc.)a. n 

Recent experiments^ with YBa2Cu30y-Nb ramp long 
Josephson junctions (LJJ) fabricated in zigzag geometry 
(if viewed from the top) clearly demonstrated that such 
a LJJ consists of alternating facets of 0, it, 0, tt . . . junc- 
tions. The presence of alternating and 7r-facets results 
in a set of 7r-discontinuities of the Josephson phase at the 
corners where 0- and 7r-facets join. These arc the points 
where the order parameter of the anisotropic high-T c su- 
perconductor changes its sign because the direction of 
the Josephson contact changes by 90°. 

In the short facet limit the I C {H) dependence of 
the zigzag LJ.J was investigated theoretically and 
experimentally^. Due to alternating and n facets the 
critical current I c has a rather small value at zero mag- 
netic field H, while the main peaks of I c are found at 
a finite field, the expression for which was derived too 
and fits well with numerical simulations and experiments 
(considering flux focusing). 

A striking property of 0-7T LJJ is the spontaneous gen- 
eration of half-integer flux quanta, further called semi- 
fiuxons (SFs), at the corners of a zigzag. .-.The presence 
of SFs was demonstrated experimentallya by scanning 
SQUID microscopy on LJJs in the long facet limit. SFs 
were also experimentally obsftriied in the so-called tri- 
crystal grain boundary LJJsffiMl In both experiments 
the samples were electrically disconnected. 

In this work we investigate numerically the behavior 
of positive and negative SFs (PSF and NSF) in a LJJ of 
zigzag geometry. The polarity of SFs is defined so that 



the PSF roughly corresponds to spontaneous magnetic 
flux +$o/2, while NSF corresponds to —<& /2. In par- 
ticular, we study the various ground states in which SFs 
of different polarity can be arranged in LJJs with many 
discontinuity points. We also investigate the effect of the 
bias current, and show that it may provoke formation of 
SFs even in the cases when the ground state has a piece- 
wise constant phase at zero bias. The rearrangement of 
the PSFs and NSFs under the effect of uniform bias cur- 
rent at zero voltage is discovered and may be observed in 
experiments similar to the one of Ref. ||. Throughout the 
paper we consider separately the cases of infinitely long 
JJ and of LJJ of finite length. Among different ground 
states we focus on antifcrromagnctically ordered chains 
of SFs (at zero bias) and on chains of unipolar SFs. We 
do not limit ourselves to only one corner but consider 
several of them, usually equidistantly spaced. 

In section [l| we introduce the model and discuss nu- 
merical aspects. Various ground states are discussed in 
section jll]. Section |lV| contains numerical results which 
show the rearrangement of SFs by means of an applied 
bias current. Section |V| concludes this work. 



II. MODEL AND NUMERICS 

The dynamics of the Josephson phase in LJJ with 
phase 7r-discontinuities is described by the following per- 
turbed sine-Gordon equationB: 

4>xx - 4>tt - sin(</>) = a<j>t - j(x) + h x (x) + 9 xx {x), (1) 

where cj>(x, t) denotes the Josephson phase across the 
junction, a = l/vT^c is the dimensionless damping 
coefficient (/3 C is the McCumber-Stewart parameteru). 
The function 9{x) describes the positions of the ir- 
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discontinuity points and can be written as 



(x) = 7T ^2 °kH(x — Xk), 



(2) 



fc=i 



where 0^ = ±1 defines the direction of the fc-th phase 
jump. The sum is over all N c corners located at x = Xk- 
7i{x) is the Heaviside step function. 

Eq. (|l|) is written using standard normalized units, i.e., 
the coordinate and the time are given in units of the 
Josephson penetration depth Aj, and the inverse plasma 
frequency ctC* 1 , respectivelyo. The subscripts x and t in 
Eq. ([j]) and below, if any, denote partial derivatives with 
respect to x and t. 

The external magnetic field enters into consideration 
through the h x term in Eq. (Q) and through the boundary 
conditions, which for overlap geometry can be written as 



where Q = tan(7r/8) = y2 — 1. The corresponding ex- 
pression for the magnetic field of a SF is 



cosh(|.T| —InQ) 



(8) 



It looks like a cusp with exponential tails. _ 

The simulations were done using StkJJ softwares fur- 
ther developed to include discontinuity points in the nu- 
merical scheme and new initial conditions with arbitrary 
distributed positive and negative SFs given by Eq. (|?|). 
For the simulations we use a linear LJJ with a damp- 
ing parameter a = 0.1, which corresponds well to the a 
values between 0.1 and 1.0 in real zigzag LJJal3. The 
applied magnetic field was zero. 



III. GROUND STATE 



cf> x (0,L) = h, 
where the field h is normalized in the usual way as 



h = 



2H 



(3) 



(4) 



Here, H c i = $ /( 7r MoAAj) is the first critical field (pene- 
tration field) for a LJJ which is, in fact, equal to the field 
in the center of the fluxon. $o is the quantum of mag- 
netic flux and A is the magnetic thickness of the junction. 
The field H in Eq. (Q) is given as H — (H • n) , where n 
is the unit vector normal to the effective cross-section of 
the junction (elementary cell). Note, that if the projec- 
tion of the external field on all facets is the same, the 
term h x (x) disappears from Eq. ([j]), and the field affects 
the LJJ only via the boundary condition (^|). 

In numerical simulations, it is quite difficult to deal 
with derivatives of (5-functions present within the 9 XX 
term in Eq. (Q). Therefore, it is convenient to present 
the total phase as a sum of two components: a mag- 
netic component /i(x) and the order-parameter related 
one 9(x) (T 



(x,t) = (i(x,t) +9(x). 



(5) 



In this case we can get the sine- Gordon equation only for 
the magnetic component of the phase /x: 

[i-xx ~ (Ht ~ sin(/i) cos(6 l ) = afn - ■y(x) + h x (x). (6) 
±i 

This is the usual perturbed sine-Gordon equation, but 
the sign in front of the sine function changes from facet 
to facet. We_use Eq. (JsJ) for numerical simulations. 

Xup^t aZ.EjO and later, but independently, Goldobin 
et ala derived an analytical expression for the SF in 
an infinitely long JJ with a single phase ^-discontinuity 
point at x = 



H(x) 



4arctan((?e :E ) ; x < 
4arctan(£e~ a: ) ; x > 



(7) 



To understand the formation of ground states in vari- 
ous configurations it is important to be able to calculate 
the energy of the system. Since the sine-Gordon Eq. (||) 
differs from the usual equation only by a factor in front 
of sin/x, the energy of the system in the general case can 
be written as 



U 



-j-oc 



fj> x + (1 — cos \x cos 9) dx. 



(9) 



Now we present the results on the observation of different 
ground states and the spontaneous formation of SFs. 



A. Infinitely long JJ 

When we speak about infinitely long JJ we mean that 
its both edges are very far from the positions of discon- 
tinuity points, fluxons, SFs and any other "interesting" 
locations. For our purposes we used L = 20 (0 < x < L) 
as a good approximation provided that SFs are at least 
5Aj apart from the edges. For some cases we repeated 
the simulations using L = 30 and L = 50, but the re- 
sults were quantitatively the same with accuracy of few 
percents. 

First, we simulated a LJJ with one discontinuity 

point at x\ = 10 and observed that the initial phase 
distribution 



(x) = 9[x) = 7tH.(x — xi), 



(10) 



is not stable and evolves into a NSF. If we take 9(x) = 
— nH(x — xi) the system generates a PSF instead of NSF. 
Note, that in principle the polarity of SF does not depend 
on the sign of the 7r-discontinuity. Both PSF and NSF 
may be pinned at a 9(x) or a —9(x) discontinuity ( [Toj ) and 
have the same energy. The energy of the SF, calculated 
by substituting Eq. (m) into Eq. (E 



U = 16- 



1 



(11) 
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while the energy of the [i = state is cx L, i.e. it diverges. 
Note that Eqs. (@) and (0) are universal for both fiuxons 
and SFs. If we take 5 = 1, Eq. converts into the sine- 
Gordon kink and Eq. (||) gives U = 8 - the rest energy of 
a fluxon. If we take Q = \[2 — 1 (as for a SF), we get the 
SF shape from Eq. (@) and U = 8 - 4^2 « 2.343 from 
Eq. (|ll). This is the rest energy of a SF. 

Second, we simulated a LJJ with two discontinuity 
points at x\ = 5 and X2 — 15. The initial phase distri- 
bution was chosen so that fi(x) — 0, i.e., 



4>{x) = 9(x) = irH.(x — xi) — irl~l(x — X2) 



(12) 



When we look at the temporal evolution of the system 
we see that such a state is meta-stable and occasionally 
degrades into two SFs of opposite polarity: a PSF and a 
NSF. The polarity of a particular SF does not depend on 
its position and on the direction of 7r jump, but depends 
on the initial perturbation which we used to disturb the 
meta-stable state. For example, when we applied a tiny 
uniform bias current 7 = 0.01, we got a NSF at x = 5 and 
a PSF at x — 15. If we apply a bias current of opposite 
polarity 7 = —0.01, we get a PSF at x = 5 and a NSF 
at x = 15. This looks natural since, for our choice of 
current direction, the Lorenz (driving) force associated 
with the positive current pushes positive (semi-) fiuxons 
to the right and negative ones to the left. Thus, the final 
configuration always corresponds to a pair of SFs pulled 
apart by the bias current. 

The states of two SFs discussed above, and naturally 
denoted as || an d T4> are stable only when the distance 
between SFs is rather large, a > a c . Otherwise, in the 
short facet limit, a < a c , the state without SFs is en- 
ergetically more favorable. The numerical simulations 

(2) 

show that the critical separation is a c ss 1.55 ± 0.05. 
The superscript stands for N c = 2, as below we also 
calculate a c for the case of more than two corners. We 
have also checked that there is no hysteresis (local en- 
ergy minimum) around a c , i.e., if we start with two SFs 
at a = a c — e (0 < e <C 1), the system relaxes to the 
/j, = state (|l2|). Inversely, if we start from the flat state 
( |l2| ) at a = a c + e, it relaxes to the f J, or JT state. 

Continuing our reasoning in a similar fashion, we con- 
clude that for odd N c in an infinite LJJ, the formation 
of at least one SF is always energetically profitable. In a 
LJJ with even N c , there is a crossover distance a c (which 
is a function of facet number, size and location) which de- 
termines the ground state. For equidistantly distributed 
corners the crossover distances a c were calculated nu- 
merically for different number of corners N c . They have 
the following values: = 1.35±0.05, = 1.15±0.05, 
a { c 8) = 1.05 ± 0.05. For odd N c the state with N c AFM- 
ordered SFs is always favored, i.e., a c = 0. In the general 
case, the situation can also be mixed, as shown in Fig. [I], 
i.e., the shorter facet A-B with a = 1 has constant <p = n, 
while the longer facet C-D with a = 2 has a PSF-NSF 
pair at the discontinuity points. 

Another case which we didn't consider up to now is 
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FIG. 1: Ground state of LJJ with 4 corners, situated in 
arbitrary way. 




FIG. 2: The local magnetic field [i x (x) of two PSFs for 
different distance a = 0.1, 1, 2 between them. 



two PSF ft (or two NSF J.J.). These states are stable 
for any separation a. For a — ► 0, the profile of two SFs 
approaches the profile of a fluxon, as shown in Fig. |^. 
The energy of the ft state changes from 2U$f ~ 4.686 
up to Up — 8 as a decreases from infinity to zero. The 
energy of the state changes from 2Usf ~ 4.686 down 
to as a decreases from infinity to zero, having U(a) w a 
for a — > 0. Although, the energy of the state is always 
higher than the one of the || state, it is still stable and 
represents a local energy minimum. 



B. Finite length LJJ 

We consider here only the case of equidistant distri- 
bution of the corners along the length of the LJJ. The 
iV c corners can be distributed in two different symmetric 
ways with the positions of the corners: 



x 



„(EF) 

'k 

(EL) 



ak, 
ak 



a 
2' 



a = 



f 1 
L 



(13) 
(14) 



where a is the facet length (the distance between cor- 
ners), and k = l...N c . In the first case, all facets 
have equal length (EF stands for "Equal Facets"), but 
for even N c the total lengths of 0- and 7r-facets are not 
equal Lj, 7^ L£. In the second case, the first and the last 
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facets are twice shorter, but the total lengths of 0- and tt- 
facets are the same for any iV c , i.e., = (EL staads 
for "Equal (total) Lengths" ) . The case of one cornertlltil 
is degenerate and belongs to both categories. 

In the analysis below we compare the flat phase state 
with the antiferromagnetically (AFM) ordered state as 
they usually represent the lowest energy solutions of 
Eq. (^). Later we also analyze unipolar states. All other 
configurations, although stable, are outside the scope of 
this paper. 

Thus, generally, the problem splits into 3 different 
cases. 

(a) Facets of equal length, even N c , Lj. ^ L^.. At 
a » 1 the ground state consists of AFM-ordered SFs, 
since it is not profitable to have long facets with <j) = ir. 
For estimation of the energy of AFM-ordered state at 
a — > 0, we note that the AFM-ordered SF chain has 

^«cos^ = (-l)'|, a^O, (15) 

where I is the facet numberEl. Then, at a — > 0, the en- 
ergy of the flat phase state U^-o = N c a, while Uafm ~ 
L = (N c + l)a, i.e., {7^=0 < Uafm- Thus the ground 
state is a flat phase state. Since we have different lim- 
iting behaviors at a — > and a> 1, there should be a 
crossover distance ai N ^ and a corresponding crossover 
LJJ length Li N "\ Our simulations for N c = 2 show 
that L ( c ] = 4.25 ± 0.05, and a { ? ] = 1.4. For N c = 4, 
L ( c 4) = 5.85±0.05, and a ( c 4) = 1.16. For N c = 6, L { c ] = 7, 
af ] = 1.0. 

(b) Facets of equal length, odd N c , L^ — L^.. |The.sim- 
plest example N c = 1 was already investigatedOEJ. At 
a ^> 1 the ground state consists of AFM-ordered SFs. 
At a — > the flat phase state and the AFM-ordered 
SF states have the same energy (in the first approx- 
imation), so itr-is reasonable to assume that for all a 
U^=o > Uafm^, and the ground state is a • • • |1| • • • 
chain. This is also confirmed by numerical simulations of 
a • • • tit ' ' ' SF chain put in LJJs of different lengths. 

(c) Twice shorter edge facets, Lf, = L|.. The ground 
state is equivalent to the previous case and confirmed by 
simulations. 

The last state which we wish to consider is tT ( or ID- 
To avoid problems with initial conditions, we have cre- 
ated a tt state in a LJJ of L = 10, and have slowly 
reduced the length of the LJJ in steps of SL = 0.1, allow- 
ing the system to relax after each small contraction. We 
have found that for each separation a there is a critical 
LJJ length L c (a), such that at L < L c two PSFs can not 
stay inside the junction. They emit a fluxon and revert 
either to TJ. state or to the flat phase state. The plot of 
L c vs. a is shown in Fig. ^. 

At a — » our two PSFs are equivalent to a fluxon. The 
pinning disappears too, as a 27r-discontinuity is equiva- 
lent to a continuous phase for our problem. A fluxon 
placed in a LJJ of any finite length is in a metastable 
state because it is attracted by two images situated be- 
hind the left and the right edges of the LJJ. Therefore, 
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FIG. 3: The dependence of minimum LJJ length at which 
the |T still exists on the distance a between SFs. 



lim L c (a) = oo, (16) 
a— »0 

Of course, in the experiment, the force pulling the fluxon 
towards the boundary is exponentially small for large L, 
and the fluxon gets pinned by non-uniformities or other 
mechanisms, provided they have a stronger effect than 
e.g. thermal fluctuations, which depin the fluxon. 

In the opposite limit a — > oo, two PSFs become in- 
dependent and L c > a, i.e., we should solve a problem 
of SF stability near the edge of a semi-infinite LJJ. This 
problem, in turn, is equivalent to the problem of stability 
of a PSF-NSF pair (the NSF appears as an image) in the 
infinite LJJ analyzed above. Therefore, the SF is stable 

(2) 

at the distance from the edge which is larger than a c /2. 
Thus, in the limit a — > oo, L c grows linearly with a, but 
is shifted by a c . The asymptote L c (a) = a + a c is 
shown in Fig. |^ and demonstrates excellent agreement of 
our arguments with numerical simulation. 

IV. REORDERING SEMIFLUXONS BY DC 
BIAS CURRENT 

An interesting question is what happens when a uni- 
form dc bias current is applied to the system. The cur- 
rent acts like a driving force which tries to push SF in a 
certain direction (just like with fluxons) which depends 
on the polarity of the SF and on the polarity of the bias 
current. On the other hand SFs are pinned by the phase 
discontinuity points and, it seems, the maximum what 
can happen is that the SF changes the shape (deforms) 
a little bit. In contrast, we found that the SF can also 
(a) spontaneously flip and generate a fluxon which starts 
moving along the LJJ under the action of a bias current, 
and (b) several PSFs and NSFs can hop simultaneously, 
if they are situated close enough, resulting in a structural 
rearrangement . 
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FIG. 4: The magnetic field fi x (x) shows the emerging SFs for 
7 = 0, 0.1, . . . , 0.7 in LJJ of L = 20 with the distance a = 1 
between corners. The pre-critical state at 7 = 0.76 is shown 
by thicker line. 

A. Infinitely long JJ 

First, we simulated a LJJ with only one discontinuity 
point. According to the results presented in the previous 
section, the ground state is a SF. Without loosing gen- 
erality let us assume that it is a PSF. When we apply 
the positive bias current from zero up to j c = 0.62 the 
PSF just changes shape. At 7 = 0.63 the PSF becomes 
unstable and splits into a positive fluxon and a NSF, so 
that the total flux is conserved. The fluxon moves away 
to the right under the action of the driving force. The 
NSF is still pinned but exposed to the same driving force 
7 > 7 C which tries to push it to the left. Correspond- 
ingly the NSF splits into a negative fluxon (anti-fluxon) 
and a PSF. The anti-fluxon moves away to the left be- 
ing driven by the bias current. Thus, after emission of a 
fluxon and an anti-fluxon, the system returns to the ini- 
tial state with one trapped PSF. Since 7 > j c , the whole 
process is repeated again, i.e., a SF of any polarity at 
7 > 7 C emits two trains of fluxons: fluxons to the right 
and anti- fluxons to the left. 

Second, we simulated the behavior of a LJJ with two 
corners at a distance a between them. When the distance 

(2") 

a < a c (we took a = 1 for simulations), the ground 
state corresponds to the flat phase profile. Nevertheless, 
if one applies a bias current, a NSF-PSF pair is formed. 
The SFs of this pair are pinned at the corresponding 
discontinuity points and the polarity of the SFs is such 
that the bias current pulls them apart. The amplitude 
(the field at the center of the PSF) of the NSF-PSF pair 
grows smoothly with the applied bias current and reaches 
the value h ps 0.75 at 7 C = 0.77, as shown in Fig. |[ This 
value of field is only twice smaller than the field in the 
center of an isolated SF, which is equal to v2 « 1.4. 
This means that even in the short face±,limit a NSF 
PSF pair can be observed experimentallylA by applying a 
uniform dc bias current through the junction. At 7 > j c , 
the NSF-PSF pair emits fluxons to the right and anti- 
fluxons to the left (for positive bias current), and, thus, 
switches to the non-zero voltage state. 

Another possibility is that the distance between two 



corners a is larger than a c (we took a = 2 for simu- 
lations). In this case the ground state can be either of 
Hi IT 7 II 7 IT; depending on the history of the system. 
Under the action of the applied bias current, the states 
|| and 11 behave very similar to a single SF state dis- 
cussed above. The most interesting thing happens when 
we apply the bias current to the || or fj, configuration, 
especially if the initial order of NSF and PSF is not natu- 
ral with respect to the direction of bias current. Suppose 
that we have the state ||. For 7 > the direction of driv- 
ing force is such that the NSF and the PSF are pushed 
towards each other, but still they can not move since they 
are pinned at the corners. Only their shape changes a lit- 
tle bit as 7 increases. At 7 = 0.08, we observe that the 
SFs flip synchronously, and exchange their location. We 

can denote such a process as T| 7 ~ —^IT- This process 
can be thought as an exchange by a virtual fluxon which 
transfers a quantum of magnetic flux from the PSF to 
the NSF: 

PSF — >NSF + Virtual Fluxon 

Virtual Fluxon + NSF — >PSF 

We use the word virtual because (a) no actual soliton is 
formed (only a quantum of magnetic flux spills over), and 
(b) this happens at the value of 7 which is much less than 
the fluxon creation current j c = 0.63 discussed above. 
The process of flipping (or hopping PSF^NSF) of SFs is 
accompanied by emission of plasma waves in both direc- 
tions, but they decay with time, so that the final state 
11 is stable. In this new state the bias current pulls the 
SFs apart and only when 7 exceeds 7 C = 0.67 the system 
switches to non-zero voltage state. If the system was ini- 
tially in the state It 7 no structural rearrangement takes 
place for 7 = . . . j c = 0.67 and the system switches to 
the resistive state at 7 = j c . Thus, by sweeping the bias 
current-hack and forth between — 7 C and +7 C one may 
observecl the structural rearrangements of SFs being in 
zero voltage state. 

Similar rearrangements we observe for the states in- 
volving many SFs. Usually we start from an array of 
AFM-ordered SFs at 7 = with the distance a — 2 be- 
tween nearest neighbors. Then we increase the bias cur- 
rent to obtain the following structural transformations. 

For 4 SFs: 

tin 7 ^ 12 mi 7 ^ 25 utt (17) 

with 7 C = 0.53 from the state ||TT- 
For 6 SFs: 

a~ 1 a~ 1 1 7=+0.13 , 7=+0.30 , , . ...... 7=+0.13 

UUU 7jL > IITITT — IIITTT 7 — ► 

iititt 7=z ^ 03 mm 7=z ^ 21 UUU (18) 

with 7c = 0.50 from the state IIITTT- I n the latter 
case the bias was first increased and then decreased, thus 
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FIG. 5: The magnetic field u x (x) profiles for different 
configurations obtained during transformation (^f|) for 7 = 
0,0.3,0.4. 



along a LJJ of finite length? For the annular LJJ the 
answer should be negative, as, due to symmetry, there 
are no left and right edges. On the other hand, for a 
linear LJJ, it seems possible. As before we distinguish 3 
different cases. 

(a) Facets of equal length, even N c . The ground state 
at 7 = can be flat for a < a^ c or consist of AFM- 
ordered SFs for a > a^ c ■ Let us start with the AFM- 
ordered state fj, at L = 6 (a=2). We observe the transi- 
tion 



ti 



7=+0.083 



IT. 



with 7 C = 0.491 from the || state. 
For N r = 4: 



(22) 



demonstrating a hysteresis in switching between various 
configurations. 

Further, for 8 SFs: 

niinn 7=±0 > 13 umm 7= - ± ^ 23 uumt 7=±0 ^ 14 
mum 7= ^ 02 nnnn 7=z °^ 20 mum (19) 

with 7 C = +0.49 from UIUTTT state. Note, that some- 
times the system switches directly to the very polarized 
state skipping several intermediate configurations. Those 
intermediate states can sometimes be obtained by sweep- 
ing the bias current back and forth, starting from various 
other states. The state UUUTT m the above example 
was obtained so. 
For 12 SFs: 

numum 7=±0 ^ 14 uumnm 7=±0 ^ 31 

UUUTTTTTT 7=± °^ 12 UUITITTTTT 7=± ^°° 

uummn 7=z °^ 05 mnnum 7=z °^ 18 

ttttUUUU 7=z ^ 3 ° TTTTTTIIIIII (20) 

with 7c = +0.47 from the state UUUTTtttt- We also 
checked that the situation is symmetric, i.e., ~/ c = —0.47 
from the state TTTTTTIIIIII- If we start from the state 
44T4ITITT4TT at 7 = obtained in J2C|), and increase 7, 
we get some new states which we have not seen before: 

IITUUTTITT 7=± °^ 12 IUUUUTTT 7= -^ 19 UUUTTTTTT 

Interesting enough, similar rearrangements are possi- 
ble even when the initial state at 7 = is flat. For 
example, for 6 SFs separated by the distance a = 1 we 
get 



7> ^°° ITITIT ^ 



HO. 37 



u-tt 

with 7c = 0.48 from the state || — |1\ see Fig. 0. 



(21) 



7=+0.119 . . . . 

UU — + IITT, 

with 7c = 0.331 from the state |UT- 
For N c = 6: 

UUU 7=± ^ 12 UTUT, 



(23) 



(24) 



with 7 C = 0.28 from the state UTITT- As we see in the 
latter case the total polarization does not take place. 
For N c = 8: 

^ , ^ , ^ , ^ 7=+0.t4 . 7=+0.23 , , , , ... . 7=+0.14 

UTITITI 7 -+ IITITITT 7 IIIITTTT 7 



IIITITTT 7 ^ 01 ITITITm 



(25) 



with 7c = 0.27 from the state UIUTTT- It is rather- 
interesting that the rearrangements happen almost at the 
same value of 7 as in the case of an infinitely long LJJ. 
We also have noticed that as a decreases towards a^" or 
as N c increases at fixed a, the current 7* at which the 
first rearrangement takes place also decreases, reaching 
7* = at a — > a^ a . In this limit the total polarization is 
almost never achieved. 

When a < a^ c the initial state is flat at 7 = 0, but 
AFM-ordered SFs emerge at 7 > 0. The amplitude of 
magnetic field at the centers of SFs grows with 7. Such 
states can not be rearranged further and usually the sys- 
tem switches to the resistive state at some 7 C . The results 
are summarized in Tab. |. The simulations also show that 
for a < the system is kind of "elastic', i.e., as soon as 
the force 7 is removed, the system, like a spring, returns 
back to the flat phase state without any hysteresis. 

(b) Facets of equal length, odd N c . The ground state 
at 7 = is an AFM-ordered . . . jj, j . . . chain, so we start 
with our usual value a = 2. The smallest odd N c for 
which it makes sense to speak about rearrangements is 
TIT: 



7=+0.t7 . . . 

TIT — > ITT, 



(26) 



B. LJJ of finite length 

Is it possible to observe similar "polarization by cur- 
rent" effects when the SFs are equidistantly distributed 



with 7c = 0.32 from the state Iff. 
For N r = 5: 



7=+0.14 7=+0.27 . . . . . 

TITIT ^ ititt — > IITTT, 



(27) 
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a 




7c 


0.5 
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0.33 


0.5 
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0.20 


0.5 


6 


0.14 


0.5 


8 


0.11 


1.0 


2 


0.36 


l.U 


4 


U.Zo 


1.0 


6 


0.17 


1.0 


8 


0.14 


1.3 


4 


0.11 


1.3 


6 


0.18 



TABLE I: The list of numerical tests which shows the 7 C 
for the switching from the ||-like states. No structural rear- 
rangements were observed (even N c ). 




coordinate, x 

FIG. 6: The magnetic field /j, x (x) profiles for different values 
of bias current 7 = 0,0.05,0.10,0.11 in LJJ of L = 6 with 
N c = 5 corners and equal facet lengths. The field of SFs is 
"pulled" to the right by the bias current. 



with 7 C = 0.28. 

For small a (we have used a = 1) the rearrangement 
of SFs does not take place, although the amplitude of 
magnetic field in the center of each SF changes so that the 
amplitude of field is higher at the right edge and lower at 
the left edge. It is like the field is globally pulled toward 
the right edge by the bias current, as can be seen in Fig. ^ 
The results of other runs for this case are summarized in 
Tab. 0. 



a 


N c 


7 C rearrangements 


0.5 


3 


0.04 


0.5 


5 


0.04 


0.5 


7 


0.04 


1.0 


3 


0.13 


1.0 


5 


0.10 


1.0 


7 


0.08 inv. left SFs 


1.5 


3 


0.17 uf^^m 


1.5 


5 


0.13 TITIT 7=± ^ 14 ITITT 


1.5 


7 


0.13 TITITIT 7= ^° 8 |TITITT 



TABLE II: The list of numerical tests which shows the 7 C for 
the switching from the J.f-like states in LJJ with uniformly 
distributed corners and equal facet lengths. No structural 
rearrangements were observed (odd N c ). 
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7 C rearrangements 


2.0 
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0.15 none 


1.0 
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0.038 none 


0.75 
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0.020 none 


0.5 
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„ „„„ , 7=+0.005 . . 

0.020 Tl ^ IT 
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0.155 none 
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0.094 none 


1.0 
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0.044 none 


u . 





0.015 none 
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0.15 none 


1.5 
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1.0 
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(J.U4U none 
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, -y — +.008 . 

0.014 till itlT 
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1.0 
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0.013 nun ^ mm 
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7 
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1.0 
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U.O 
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7 
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1.0 
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0. none 
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0.012 umm 7=± ^ 09 mum 



TABLE III: The list of numerical tests which shows the value 
of 7c for the switching from the ||-like states in a LJJ with 
uniformly distributed corners and twice shorter edge facets. 
No structural rearrangements were observed. 



(c) Twice shorter edge facets. The ground state at 
7 = consists of AFM-ordered SFs. Numerical results 
show that no rearrangements were observed for iV c = 
8, . . . , 18 in a LJJ of L = 30 (a from 3.75 to 1.67). For 
smaller a rearrangements were observed in few cases at 



0, see Tab. [II 



Such a difficulty of rearrangement 
can actually be explained. If we remove the LJJ edges 
and introduce SF images, our problem reduces to the 
problem of rearranging an infinite chain of AFM-ordered 
SFs. Similar to the case of an annular LJJ, there is no left 
and right ends here and all points are equivalent, so the 
rearrangement is never initiated. Synchronous hopping 
of all SFs would require infinite energy. Of course, in 
an experiment the facets are never perfectly equal, so 
this may result in rearrangements of SFs. In simulations 
rearrangements may take place due to the finite Aa; of 
the numerical scheme or the finite step in 7. At a — > 
the rearrangement, if any, happens at 7 — > 0. 

The hopping of SFs under the action of bias current 
results in a possibility to observe half-integer Zeso-Field 
Steps (ZFS) in a LJJ with moderate facet sizeEa. If we 
take a LJJ of length L = 2 with a phase discontinuity 
point and a PSF in the middle, the interaction of the 



8 



PSF with the boundaries at zero applied magnetic field 
can be treated by introduction of two NSF (images) at 
the distance 1 from the left and right edge of the LJJ. 
In this case the boundary can be removed and we have 
a pinned NSF-PSF-NSF system with a = 2 in an infinite 
LJJ. Under the action of a bias current 7 > 7*, the PSF 
exchanges its position with one of the images (NSF) so 
that we get a NSF in the center of LJJ and one quantum 
of magnetic flux (virtual fluxon) transferred through the 
edge. Now the system is equivalent to PSF-NSF-PSF 
with a = 2 in an infinite LJJ and the whole process is re- 
peated again, transferring another virtual fluxon through 
the edge of the LJJ. At this point the state of the system 
is exactly equal to the initial state and the whole period 
is repeated again. By analogy with ZFS, the PSF is also 
"reflected" as an NSF after interaction with the edge, but 
this process is not continuous as in the case of fluxons, 
but somewhat discrete — a PSF hops out and a NSF 
hops in. Since the total flux transferred per "reflection" 
is 2w (twice smaller than for conventional ZFSs), and the 
velocity of the virtual fluxon presumably does not exceed 
the Swihart velocity one should observe half-integer ZFS 
in the experiment. 

V. CONCLUSIONS 

We have discussed various configurations of semi- 
fluxons (SFs) in long Josephson junctions (LJJs) with 
7r-discontinuities (zigzag corners) and corresponding 
ground states. We have shown that (a) in an infinite 
LJJ the ground state depends on the separation between 
the corners (facet size), and may be either a flat phase 
state, or an array of antifcrromagnetically ordered SFs, 



or an array of unipolar SFs. If the number of corners is 
odd, at least one SF should be present, (b) In a LJJ of 
finite length the ground state of a LJJ with N c uniformly 
distributed corners can be flat or consist of SFs depend- 
ing not only on the facet size a, but also on the question 
whether iV c is an odd or even number and on the way 
to distribute the corners, (c) In the case of a flat phase 
ground state, one may still compel SFs to emerge by 
applying a uniform dc bias current, (d) The rearrange- 
ments of SFs by bias current, .discussed in this paper, 
can be observed experimentalljO by being in the zero- 
voltage state and changing the current in the range from 
—7c to +7 C . One may observe the whole variety of states 
discussed here, if the number of facets and their sizes are 
properly selected, (e) Although the SFs are pinned, they 
may hop from one discontinuity point to the next one, 
provided the distance between these points is about 2Aj 
or less. This "moving by hopping" results in a possibility 
to observe some phenomena known from fluxon dynam- 
ics, e.g., zero-field steps, but with the first step situated 
at half of the voltage of the usual first zero-field step. 

The most numerical observations reported here can 
be directly checked in a SQUID microscopy experi- 
ment. The results presented here provide new insights 
into the statics and dynamics of SFs in LJJs with n- 
discontinuitics. We also hope that this work will stimu- 
late further research and generate ideas for novel devices 
and their application in superconducting electronics. 
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